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Displacement field and elastic constants in non-ideal crystals 
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In this work a periodic crystal with point defects is described in the framework of linear response 
theory for broken symmetry states using correlation functions and Zwanzig-Mori equations. The 
main results are microscopic expressions for the elastic constants and for the coarse-grained density, 
point-defect density, and displacement field, which are valid in real crystals, where vacancies and 
interstitials are present. The coarse-grained density field differs from the small wave vector limit of 
the microscopic density. In the long wavelength limit, we recover the phenomenological description 
of elasticity theory including the defect density. 

PACS numbers: 62.20.D-, 46.05. -fb, 61.72.jd, 61.72.jj, 63.20.-e 



I. INTRODUCTION 

The theory of elasticity of solids started with Hooke 
in 1678, when he formulated the linear relation between 
stress and straii*^. The atomistic picture of matter con- 
tributed a quantitative microscopic understanding of the 
mechanical properties of ideal crystals based on the parti- 
cle potentials^. Yet, the restriction to ideal crystals con- 
taining no point defects needs to be stressed. Noncqui- 
librium thermodynamics achieved a phenomenological 
description of the long wavelength and low frequency 
excitations^i^. Martin, Parodi, and Pershan showed 
that the spontaneous breaking of continuous transla- 
tional symmetry leads to eight hydrodynamic modes, one 
of which corresponds to defect diffusior^^. Point defects, 
like vacancies and interstitials, arc present in any equilib- 
rium crystal and a complete microscopic theory of crys- 
tal dynamics needs to include them. Interestingly, such a 
complete microscopic theory of real crystals was lacking 
and is developed in this contribution in the framework of 
linear response and correlation functions theory of broken 
symmetry phases. 

Crystals exhibit long-range translational order and 
possess low-frequency Goldstonc modes, e.g. transverse 
sound waves, which try to restore the broken symme- 
try. In the familiar microscopic description of ideal crys- 
tals, the long-range order is incorporated at the start 
by assuming that the equilibrium positions of the parti- 
cles are arranged in a perfect lattice. A one-to-one map- 
ping follows between the particle i and its lattice position 
R* = (r'(t)). The deviation between the actual position 
v^{t) and the lattice position R* is called displacement 
vector u*(t) 

u'(t) = r'(i) - R' (ideal crystal). 

The (symmetrized) gradient tensor of the displacement 
vector field is connected to the strain tensor, which plays 
the central role in the theory of elasticity. Yet, the appli- 
cability of the displacement vector is, due to the one-to- 
one mapping, restricted to perfect crystals, because an 
interstitial corresponds to a particle without lattice site 
(Fig.[T^), and an vacancy to a lattice site without particle 
(Fig.[T}D). Moreover, because defects are mobile, any 'im- 




FIG. 1: Schematic two dimensional drawing of a crystal con- 
taining (a) an interstitial and (b) a vacancy with threefold 
symmetry^. 



proved' mapping would yield displacement vectors that 
can become arbitrarily large with time. Linear elasticity, 
considering small strain fields, thus would intrinsically 
be restricted to short times, contradicting/invalidating 
its application to low frequency vibrations. Thus the 
need arises to define the displacement field microscopi- 
cally without the recourse to a perfect lattice of equilib- 
rium sites R*. 

In an ideal crystal where the one-to-one mapping of 
particles to lattice positions holds, a density change is 
given by the divergence of the displacement field^i^ 

5n{Y,t) = — tt-qV • u(r, t) (ideal crystal), 

with average density no = N /V and obvious definition of 
u(r, t) = u*(i)(5(r — R*); note that it is a periodic lat- 
tice function in this case. Above relation holds because 
the density can only fluctuate by particles moving around 
their lattice sites. In an ideal crystal, density thus is not 
an independent degree of freedom, and description of the 
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displacement field suffices. In a real crystal containing 
point defects, translational symmetry is still broken and 
long range order exists, but the motion of defects de- 
couples density fluctuations and the divergence of the 
displacement field. Density not only changes because of 
a deformation of the lattice, given by —V • u{r,t), but 
also by motion of additional/missing particles from one 
lattice cell to another. The density can be decomposed 
in two parts: 

6n{r,t) = ~noV ■u(r,t) ~ Sc{r,t) . (1) 

This definition of the defect density c(r, t) is positive for 
vacancies, and negative for interstitials, i.e. the conven- 
tional sign favors the interpretation as a vacancy density. 
Comparison with the discussion in Q for vacancies as the 
only point defects shows, that the thus defined variation 
in defect density is given in terms of number of vacancies 
A^^ and number of lattice sites Nis 

Thus, the magnitude of the variation of vacancy density 
increases, if there are more vacancies, and decreases, if 
there are more lattice sites. 

Based on the relation ^ alone, the hydrodynamic pre- 
dictions by Martin et al^ can be recovered. Yet, there ex- 
ists no microscopic particle based theory, which provides 
the definitions of the displacement and defect density 
fields, and recovers Eq. ([T]) from first principles. We will 
present these definitions and derive the equations of mo- 
tion for the fields, which reduce to the continuum descrip- 
tion in the hydrodynamic limit. We will follow, within 
linear response theory, the accepted route to symmetry 
broken states by considering conserved and symmetry 
restoring fields — based on an application of Bogoli- 
ubov's inequality — , followed by Zwanzig-Mori equations 
as pioneered by Kadanoff and Martin^, and Forster— , and 
by taking the hydrodynamic limit at the end. 

In an important contribution, Szamel and Ernsli^ sug- 
gested the definition of the displacement field that we 
will find, which only uses density measurements with- 
out recourse to an underlying lattice. Importantly, the 
new expression for u(r,t) can thus be used for both 
ideal and real crystals, either by simulation or, experi- 
mentally, by optical techniques in e.g. colloidal crystals. 
Because we have a systematic way to discover the hy- 
drodynamic fields and their equations, we can correct 
the work by Szamel and Ernst and achieve consistency 
with the phcnomcnological description, which these au- 
thors could notiSiii. Our approach uses density func- 
tional theory (DFT), to describe the equilibrium correla- 
tions in a crystal, and thus superficially bears similarities 
to earlier works using approximate DFTsi^"— . In con- 
trast to these previous works, we use exact DFT relations 
to simplify our expressions, and do not approximate the 
free energy functional, nor start from parametrizations 
of density fluctuations^; see Kirkpatrick et al^ for a 



discussion of these approximate theories, and computer 
simulationaiSi^i22ri2i for possible problems arising con- 
cerning the elastic constants^^. 

The paper is organized as follows: Section II de- 
rives the Zwanzig-Mori equations for (classical) crys- 
talline solids, where translational symmetry is sponta- 
neously broken and long range order exists. For simplic- 
ity, the set of conserved variables is restricted to density 
and momentum, neglecting energy. This restricts us to 
an isothermal approximation. Again for simplicity, mem- 
ory kernels are neglected, restricting us to a dissipation- 
less theory. Because the complete (infinite-dimensional) 
set of symmetry restoring variables, derived from Bogoli- 
ubov's inequality, is considered, a systematic approach to 
the dynamics of crystals is achieved; because we use the 
fiuctuation dissipation theorem, a theory linearized close 
to equilibrium is obtained. Section III identifies the con- 
ventional fields used for describing the dynamics of crys- 
tals, especially the displacement and defect density field. 
Their equations of motion within the first Brillouin zone 
are derived. Section IV uses symmetry considerations 
within density functional theory, to derive the properties 
of the coefficients entering the equations of motion, and 
Sect. V discusses the results. First, the phenomenologi- 
cal equations of elasticity theory are recovered, and the 
elastic constants identified; we obtain their microscopic 
expressions in terms of the direct correlation functions 
of the crystal. Then the displacement and defect density 
field arc discussed. Section VI ends the main text with 
short conclusions, and Appendix A shows consistency of 
the conventional but simplified Zwanzig-Mori equations 
of a crystal to our results. 

II. GENERAL THEORY 

A. Microscopic Model and Microscopically Defined 
Hydrodynamic Variables 

We consider a volume V containing TV identical spher- 
ical particles at number density uq = N/V. The motion 
of the particles with identical mass m is described by 
a (classical) Liouville operator £, which includes kinetic 
and (internal) potential energies. 

For the derivation of hydrodynamic equations, the con- 
served quantities need to be considered. Starting with 
particle number, the (fluctuating) microscopic density is 
a sum over all particles i 

N 

p(r,t) = ^<5(r-r^(t)). (2) 

1=1 

Temperature T and density no are chosen such, that the 
crystalline state gives the lowest free energy and transla- 
tional invariance is spontaneously broken. Long-ranged 
order exists and the average density varies periodically 

n(r) = (p(r,t))=^nge's-r^ (3) 

g 
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where the order parameters rig are the Bragg-peak am- 
pHtudes at the positions of the reciprocal lattice vectors 
g, which are defined by 

g • L = 2™ , (4) 

where n is an integer, and L the set of discrete transla- 
tional symmetry operations in real space. This means 

7i(r) = n(r + L) V L . (5) 

An ensemble of identical crystals, which are just dis- 
placed in their center of mass or overall orientation, yield 
vanishing order parameters. To specify the broken sym- 
metry state, it is thus necessary to fix the six degrees of 
freedom of a rigid hody^. Conceptually one describes 
the system in a frame of the center of mass and orien- 
tation, or confines the crystal with the help of external 
potentials. An example of such a potential is an exter- 
nal wall, which in thermal equilibrium would need to 
be placed such that Nis crystal lattice sites fit into the 
volume without externally applied macroscopic strain or 
stress; here Nig differs from TV because of point defects 
like vacancies and interstitials. The (canonical) ensem- 
ble, used to define the averages (...) = J p. . .dT, and 
the corresponding C is henceforth restricted to contain 
such a device which fixes the degrees of freedom of a 
rigid body. Because the internal fluctuations are not af- 
fected in the thermodynamic limit, our results will de- 
pend on the canonical set of thermodynamic variables 
(temperature T, number density rtg, volume V), and the 
order parameters rig . Because they take their equilibrium 
(non-strained) values, our (later) use of the fluctuation- 
dissipation theorem restricts us to obtain the linear equa- 
tions of elasticity, linearized around the equilibrium at 
vanishing displacement field, (u) = 0. 

The standard Fourier transformation in space is used, 
where d is the spatial dimension, and it will be stated ex- 
plicitly if a specific spatial dimension is considered, which 
usually will be three dimensional space. 

r ^ 

p(k, t)= d''re-'^-^p{r, t)^J2 e"'*""""^*^ , (6a) 

Here the reciprocal vector k is unrestricted. The lattice 
symmetry also leads to periodicity in reciprocal space, 
which can be considered to be composed of periodically 
arranged Brillouin zones. If one restricts the recipro- 
cal vector q to the first Brillouin zone, then the Fourier 
transformation of the density can be unambiguously de- 
composed into a reciprocal lattice vector g and q 

p(g + q, = / rf'Ve-^(s+q) r^(j.^ ^) _ (7) 
The Fourier-back transformation simply becomes: 



This splitting of the Fourier coefficients of the density 
is useful as for the hydrodynamic description one is in- 
terested in the long-wavelength fluctuations, i.e. g — > 0, 
close to all positions g of the order parameters n^. Using 
the Fourier-transformed density the rig are identified as 

1 1 ^ 

i 

The second conserved quantity, to be considered is 
momentum. For the momentum density ja(r,i), which 
straightforwardly is given by 

N 

Uv,t)^Y.^J{v^v\t)) , (10a) 

i 

j^(k,t) = J d^e-^'^-^Uv^t) = Epie-'''--'« , (10b) 

i 

the distinction between k and q is not necessary. (Greek 
indices are used for spatial components, whereas latin 
ones denote particles.) 

The conservation of particle density is expressed via 
(use of Einstein's sum convention is implied) 

m9tp(k,i) + i/c„j„(k,t) = 0, (11) 

which follows from the microscopic definitions Eqs. © 
and (jlOp . The conservation of momentum density is 
stated through the divergence of the stress tensor. 

9tj„(k,t) - %cr„/3(k) = . (12) 

A microscopic definition of the stress tensor can be found, 
for example, in As third conserved field, the energy 
density should be considered. For simplicity however, we 
neglect the coupling of energy fluctuations to the mechan- 
ical fluctuations, restricting our results to an isothermal 
approximation. Extensions, incorporating energy fluctu- 
ations are straightforward, in principle. 



B. Bogoliubov Argument 

In a state with spontaneously broken symmetry, ad- 
ditional variables besides the conserved quantities need 
to be considered for deriving the continuum mechanics 
equations. This by now classical route to 'generalized 
hydrodynamic or elasticity theory' — in contrast to 'hy- 
drodynamic theory without broken symmetry' — builds 
on the Bogoliubov inequality to identify variables with 
long-ranged equilibrium correlations. For crystals this 
variant of Schwarz's inequality has been formulated by 
Wagner—, 

/|c .„^„^|2^ . |(j:(k)/:<5p(g + q))|^ 
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where use is made of the hermitian property of {{CA)* B) , 
and 5p denotes the density fluctuations from the equihb- 
rium density 

5p{v,t)^ p{Y,t)-n{T) . (14) 

The correlation functions required for the Bogoliubov in- 
equality are considered in the following. Most of them 
will also be useful for elements of the so called frequency 
matrix. First the classical equipartition theorem states 
for the correlation of the different spatial components of 
the momenta of the particles 

{PaP'li) = mkBTSijSap ■ (15) 

With this 

(<5j: (k, t)6ji3 (k, t)) = mnoVkBTdc.13 ■ (16) 

The standard properties of the Liouville operator—, here 
C = —idt, yield 

m 

(17) 

Using Eq. ^ the numerator of the Bogoliubov inequality 
becomes 

{j:ik)CSpig + q)) =-^^^ 0:(k)jMg + q)) 

= -(g + 9)afcBT^^e-'(s+q-k).r'^ 

= -{9 + q)akBTVn^+q-k ■ (18) 

Thus k has to differ from q in the first Brillouin zone by 
a reciprocal lattice vector g' in order to give a finite order 
parameter ng„g'. The denominator is expressed with the 
conservation of momentum density, Eq. (|12p . 

(l'CjQ(k)P) = kp{(Tap(Taj)k^ , (19) 

or, neglecting directional dependence, ([T^ — Rk^ with R 
the correlation of the stress tensors. 
So finally 

(lM.+q)l->> '^^'"^"^■'g;^^--|'^^ (20) 

As 7ig+q-k 7^ only ifg + q — k = G, with G a recipro- 
cal lattice vector, the Bogoliubov inequality becomes 

(lMG + k)p>><£±i£<^^g>>^iPX!„*-^ 

(21) 

Note that 6p{G + k) is well defined for fixed G, and in 
the hydrodynamic limit k can be replaced with q, 
which lies in the first Brillouin zone. 

Importantly the Bogoliubov inequality is an argument 
for all G 7^ 0, but not for G — 0, as in this case the 



right hand side of ((2T|) is proportional to This com- 
ponent of the density is just conserved, it does not refiect 
the broken symmetry. Otherwise, at all finite reciprocal 
lattice vectors, the (expected) Bragg peak, which arises 
from the coherent scattering and is infinitely sharp in 
the present treatment because of the long-ranged order, 
sits on top of a diverging (diffuse) background. Bogoli- 
ubov's inequality only gives a bound for the divergence 
for A: — > 0. In Sect. IIVB II we will apply relations from 
density functional theory to prove the vanishing with k^ 
of re-summed elements of the inverse of the density cor- 
relation functions, which corresponds to the equality sign 
in Eq. 

Following the standard reasoning to derive generalized 
elasticity theories, the 'symmetry restoring' fluctuations 
need to be included in the set of slow variables and 
lead to Goldstone modest. Equation (PT|) shows that 
this requires to include the density fluctuations close to 
all reciprocal lattice vectors G; note that we will apply 
Eq. (PT|) with the trivial notational change, G replaced 
by g, in the following. 

C. Zwanzig-Mori Equations of Motion 

Whenever a set of slow variables {Ai{t)} is selected, 
the Zwanzig-Mori formalism yields their linear equations 
of motion^ '^^i^^ . Neglecting dissipation, i.e. memory ker- 
nels, the 'reversible' equations of motion for small devi- 
ations (SAilt))^" are given in terms of the equilibrium 
frequency matrix flik 

5t(M,(t))'^-ia*,(Mfc(t))''- (22) 
- i[{SA*CSA,){6A*SAk)-^y{6Ak{t)y'- . 

Here following Onsager and the fluctuation dissipation 
theorem, the deviations {SAi{t)y'^ of the spccifled vari- 
ables from their equilibrium values are within linear re- 
sponse connected to correlation functions evaluated in 
the unperturbed system. The averages and the Liouville 
operator C defining the frequency matrix in Eq. (|22|) thus 
belong to the canonical ensemble introduced in Sect. Ill Al 
The equations of continuum mechanics can be derived 
from Eq. ((22)) by choosing as slow variables the set of 
conserved and broken-symmetry restoring densities, and 
then analyzing the limit of small wave vectors, q . 
Based on Bogoliubov's inequality Eq. (PT|). the set of vari- 
ables comprises the d components of the conserved mo- 
mentum density 6ja{c{, t), and the Fourier components of 
the density fluctuations Sp{g + q, t) close to the Bragg- 
peak positions; to uniquely specify the latter, let us recall 
that the wave vector q is restricted to lie in the first Bril- 
louin zone. To clarify the notation in the following, we 
abbreviate 

dn^{q,t)^{Spig + q,t)y\ (23a) 
<5ja(q,t) = (<5ja(q,i))'' ■ (23b) 
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Zwanzig- Mori's equations (|22|) will for this choice of vari- 
ables in the limit of small wave vector lead to the dissi- 
pationless, isothermal (i.e. neglecting coupling to heat 
flow), and linearized equations of crystal elasticity. 

Most of the elements of the frequency matrix ilik have 
been derived in the previous section, namely in Eqs. (|16|) 
and (fT5|) ; note that the latter will be used in the following 
for wave vectors k = q in the first Brillouin zone only. 
Many matrix elements vanish because of symmetry^. In 
the problem at hand the most useful symmetry in this 
respect is invariance under time reversal, as the dynam- 
ical variables have a definite parity (even for the Fourier 
components of the density and odd for the momentum 
density), as well as the Liouvillian C (odd). 

The only non- vanishing matrix element still missing is 
the inverse Jgg' of the density correlation matrix, which 
is defined by 

VkBTd^s" = Y.(^P*{g + CDS Pis' + q)) Jg'g"(q) • (24) 
g' 



Its properties will be discussed in Sect. IIV Al 

Thus the dissipationless and isothermal Zwanzig-Mori 
equations of motion of a crystal are 



'-^'^-^^'^^^^l J^.-(q)%(q)) j^^^^^'^^ 

771710 



(25a) 



g'.g 



{Sp*{g' + q)Sp{g + q)) 



= -i^ig' + q)anl,J*,g{q) (57ig(q,t) 



(25b) 



Although formally exact, these equations still need a lot 
of interpretation. To begin with, there are oo + 3 of 
them in three dimensions. Naturally the question arises 
how the set of {(57ig(q, t),6ja{(l, t)}, in the limit of small 
wave vector q, reduces to the seven conventional ones 
of elasticity theory {Sn{c[,t),Sua{q,t),Sja{q_,t)}, if the 
coarse-grained density Sn(q,t) (or instead the vacancy 
density 6c{q,t)) and the displacement field 6ua{c[,t) are 
used, as in the case of phenomenological theory^i^ii (see 
Sec. IV Bl for a summary of phenomenology). In terms 
of the frequency matrix this corresponds to solving the 
eigenvalue problem, thus showing that this matrix has 
seven eigenvalues which become arbitrarily small in the 
limit g — >■ 0. These eigenvalues are the ones of classi- 
cal elasticity theory, and their corresponding eigenvec- 
tors are the variables of the continuum approach derived 
within our microscopic theory. 



III. RELATION TO CLASSICAL ELASTICITY 

The Zwanzig-Mori equations of motion (|25|) can be 
written with the frequency matrix in a compact nota- 



tion, in order to analyze the hydrodynamic limit. Before 
doing this in general, the wave equation, which contains 
the constants of elasticity or sound velocities, can be read 
off immediately. 



A. Wave Equation 



Taking a time derivative of Eq. (j25bp , and combining 
it with Eq. (|25ap . leads to a closed equation of motion 
for the momentum density 

d^SMq,t) = 

— YXg' + 9)a<, J* (q)7ig(.g + q)i3Sjp{q,t) 



771710 

g'.g 



mno 



^apiq) Sjp{q,t) . 



(26) 



This equation can take the required form of the wave 
equation, if the d x d-dimensional matrix Kap (q) vanishes 
quadratically with wave vector going to zero, AQ,^(q) = 
0{q^) for g — > 0. This property and the relation with 
the elastic constants, which obey the Voigt symmetry 
in their indices, is the subject of chapter IIVI Strictly 
speaking, only then the term wave equation is justified. 
From Eq. ([25]) one reads off 



Aa^i(q) = Y^{9' + g)ang'>/g'g(q)»^g(,9 + q)fi ■ (27) 
g'.g 



The remarkable feature of this equation, however, is 
that it is exact and holds for wave vector q throughout 
the first Brillouin zone. It is independent of the yet 
to be found relation between (57ig(q, i) with the dis- 
placement field 5ua{q^t) and the defect density 5c{q,t). 
The only input for this relation are the exact matrix 
elements of (<5j* (q)%(q)), {jl{q) C5 p{g + q)) , and 
(,5p*(g + q)<5/i(g'+q)). 



B. Displacement Field and Defect Density 

The Zwanzig-Mori equations of motion (|25p of the set 
of conserved and Goldstone modes couple density fluctu- 
ations with modulation given by (almost) the reciprocal 
lattice vectors g. To bring out the contributions from the 
various g, consider 5fi{q, t) as infinite-dimensional (col- 
umn) vector, whose components (57ig(q, t) are indexed by 
the g (ordered in some fixed but arbitrary fashion). Let 
Va,% = i{g + <l)ang be an element of a constant infinite- 
dimensional vector Va (index a continues to label the 

spatial coordinate), and let J~ J^^i be a correspond- 
ing oo X cxD-dimensional matrix. Thus the Zwanzig-Mori 
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equations (|25l) take the form 



dtSn{q, t) 



1 



mno 



(28a) 



dtS Jain, t) = J Sn{(i,t) , (28b) 

with = —i{g + q)ang an infinite-dimensional (row) 

vector, and J 5n a shorthand notation for J*^,5ng/ , 
etc. 

The prequel of the wave equation, Eq. (j26p , of the mo- 
mentum density immediately is recovered and takes the 
form 

2 1 

dtSja{(l,t) = J vpdjp{ci,t) = A„^(5j^(q,i). 

mno muQ 

(29) 

In a perfect crystal, density fluctuations result from 
the divergence of the displacement field^*^ 

Sn{q,t) ~ ~inoqadua{q,t) , (30) 

motivating the following defining relation 

(5n(q, t) ~ —VaSua{c[,t) , (31) 

for the (Fourier transformed) displacement field 
6ua{ci,t). Equation (|25a[) then becomes 



Va \ dtSua{q,t) 



(32) 



which states for an ideal crystal, as expected, that the 
time derivative of the displacement is the velocity field, 
i.e. the momentum density field divided by the mass den- 
sity 



dtdua{q,t) 



1 



mno 



(33) 



Generalizing this consideration in a real crystal, the 
(Fourier transformed) defect density field Sc{q, t) can be 
defined by the difference between the density fiuctuations 
and the divergence of the displacement field 



Sn{q, t) = -Va6ua (q, t) + v'6c{q, t) , 



(34a) 



with v' some yet unknown constant vector with com- 
ponents Wg. For convenience this can be rewritten with 

Va = Va{q = 0), whose components = icjang are the 
Fourier components of the gradient of the equilibrium 
density. 

(Sn(q, t) = -v^6ua (q, t) - iqafiSua (q, i) + v'Sc{q, t) , 

(34b) 

with n the vector with components Tig. It is an ansatz 
that the four (i.e. d -I- 1) dynamical variables introduced, 
namely the ^^^(q, t) and Sc{q, t), together with the three 
(rf) familiar components Sja (q, t) , solve the infinite set of 
equations (|25p . whose uniqueness we can not prove. Yet, 



the proof that this ansatz solves Eqs. ((25)) . is straightfor- 
ward and leads to relations expected from phenomenol- 
ogy (see Sec. lVB|) . Entering Eq. ((Ma|) into Eq. (P5a|) and 

decomposing = + iqan, one arrives at 



Va \dtSua{q,t) —Sja{q,t)\ 

L mno J 



mno 



-qaSja{q,t) + iqadt6ua{q,t) \ - v'dtSc{q,t) . 



If one continues to require that the time derivative of the 
displacement is the velocity field, Eq. (|33|) , then the first 
bracket vanishes as before. As for the general case with 
dissipation dtSc{q, t) one needs a further relation for 
the vanishing of the right hand side for all q. This is 
achieved by taking the unknown constant vector v' = 
—n/no- One is then able to define the coarse-grained 
density variation by the expected relation ([1]) 

6n{q, t) = -inoqaSuaiq, t) - (5c(q, t) , (35) 

which states that density fiuctuations are composed of 
the divergence of the displacement field and defect den- 
sity fluctuations. Consequently, the original Eq. (|25ap is 
solved for all q by the conservation law of mass or particle 
number, which causes also the second bracket to vanish 



dt8ua{q-,t) 



mno 



-5ia{q,t) 



=0 



-n —dt5n{q,t) H —qaSja{q,t) 

no mno 



Equation (|25ap thus is solved by the ansatz for all q in 
the first Brillouin zone. 

Turning to the second Zwanzig-Mori equation (j25bp . 
the ansatz Eq. (|34ap transforms it into 



, <->* , <->* n 

dtSja{q,t) = -v^ J Vf!5up{q,t) ~ vl J — 5c{q,t) 

no 

= ~kap5up{q,t) - Va(5c(q,t) , (36) 
with constant d-dimensional vector Va given by 

Va{q) =vl J — = — - y^ig' + q)an*,J*, {q)ng . 
no no ^-^ 
s',s 

(37) 

This equation is consistent with the wave equation for 
the momentum density, when a time derivative is taken, 
Eq. (|33|) is used, and the defect density is constant in 
time 

- Vadt Sc{q, t) ^ d^6ja{q, t) + -^KafsSjpi^, t)^0. 

mno 

(38) 
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Alternatively, Eq. (|36)) leads to the wave equation for 
the displacement field when the time derivative of the 
displacement is again identified as velocity field 

dtSuaiq, t) = 3_A„^(q)fe^(q, t) '^—Va{q)5c{c{, t). 

mnQ mriQ 

(39) 

Importantly, the matrix of elastic coefficients A^^ (q) in 
this prcquel of the wave equation is identical to the one 
derived for the momentum density, Eq. (|26p . Equation 
([55)1 can thus reduce to the expected result from classi- 
cal elasticity theory in the limit of small wave vector, if 
Va{q) = 0{q) can be shown for g — )■ 0. This will be dis- 
cussed together with the properties of A^^ in Sect. lIVBT] 
Then we will be able to conclude that the Zwanzig-Mori 
equations of a crystal (^5)) are solved by seven {2d + 1) 
coarse-grained fields which are the momentum density 
Sja (q, t) , and the displacement Sua (q, t) and defect den- 
sity (5c(q, field introduced in Eq. (|34a|) . We will also 
be able to conclude that the coarse grained fields, whose 
equations of motion were just determined in this sec- 
tion for q in the first Brillouin zone, obey the equations 
of motion known from phenomenological elasticity the- 
ory in the limit of g — > 0. In support of this, Eq. ([5^ 
finds that the prefactors in front of 6up{q, t) and (5c(q, t) 
are connected. Their q 0-limits reduce to thermo- 
dynamic derivatives, which, from equilibrium thermody- 
namics, have to satisfy Maxwell relations^. While the 
relations are not sufficient to express one coefficient in 
turn of the other, the connection between Aq,^ and Va in 
Eq. p9p closely mirrors the thermodynamic one expected 
in classical elasticity theory; see Sect. IV Bl 

The question remains how 5ua{q^^ t) and (5c(q, i), given 
implicitly in Eq. p4ap , can be obtained directly in terms 
of the density fluctuations (5ng(q, i). Fortunately, this 
can be achieved by projecting (5n(q, t) onto the two vec- 
tors and n, as they are orthogonal 

« (^a)^^ = El"gl'5a=0, (40) 
g 

because of symmetry. Projecting the ansatz for the den- 
sity fluctuations in Eq. p4ap onto gives an explicit 
formula for the displacement field 5ua. (q, t) in terms of 
the 5ng(q, t) 

5ua{(\, t) = lAf^i^ E 71* gp (5ng(q, t) , (41) 
g 

with Afap ~ X]g V'^^QaQfi- A second summation over 
(5?'ig(q, t) obtained from projecting Eq. p4ap on n yields 
the hydrodynamic variation of the coarse-grained density 

Mq,i) = ^E"g^"g(q'^)' (42) 

^ g 

where A/q = X]g I'^gP- With ([55]) . the relation between 
the variation of coarse-grained density, the lattice den- 



sity, and the defect density, we get 

(5c(q, t) = -no X! "g [jf ^ ■'^apladis) ^%(q, t) ■ (43) 
g ° 

IV. SYMMETRY AND INVARIANCE 

This chapter completes the derivation of Zwanzig-Mori 
equations of motion in terms of the Fourier components of 
the density by showing that the proper characteristics are 
recovered in the hydrodynamic limit. Properties of the 
(inverse) density correlation matrix lead, due to transla- 
tional invariance, to the correct q dependence and, due 
to rotational invariance, to symmetries of the constants 
of elasticity. 

A. Symmetries of Density Correlation Functions 

The symmetry property in Eq.Q of the average den- 
sity of a crystal is very familiar. Also important is the 
symmetry property of the equilibrium two-point correla- 
tion functional. For example the correlation of the den- 
sity fluctuations 5p{r, t) — p{r, t) ~ n{r) is also periodic 

{6p{ri)5p{r2))^{6p{ri+L)6p{r2+L)) VL. (44) 

This results in a periodic center of mass variable R = 
'"^^''^ and a Fourier coefficient, which depends on the 
difference Ar = ri — r2 

(<5p(ri)<Sp(r2)) = Ee's '^fc(2)(ri - r^) . (45) 
g 

As the density is a real quantity and the correlation func- 
tion is symmetric with respect to interchange of variables, 
the Srig '{rx — r2) obey the following two equations 

fc(.2)(ri - ra) = (5nL^^(ri - ra) = dn^^Hr2 - ri) . (46) 
Rewriting Eq. (|45|) , one realizes that the Fourier transfor- 

(2) 

mation of the Fourier coefficient Srig ( Ar) with respect 
to the difference coordinate Ar can be understood as a 
generalized structure factor S'g(k) 

5g(k) = 1 y d-Viy d'Va (Jp(ri)(5p(r2)) e-'^'^ e-''^'^'- 

(47a) 

= y'd'^Ar (5n(.2)(Ar)e-*-^'" (47b) 
= dn^^\k) = ^{Sp{g/2 + k)5p{g/2 - k)) . 

(47c) 

The generalized structure factor S'g(k) is, due to the sym- 
metry of the crystal, a density fluctuation function eval- 
uated with a combination of a reciprocal lattice vector g 
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and a reciprocal vector k. Bogoliubov's inequality shows 
that So{k) diverges quadratically at aU reciprocal lattice 
vectors k — )• g 7^ 0, which however is not enough informa- 
tion to simplify completely the expressions for the elastic 
coefBcients that we derived in the previous Sect. IIIII It 
remains to study the complete matrix of inverse density 
correlations, Jgg/(q) defined in Eq. (P^ . which we under- 
take now using density functional theory^"— together 
with the symmetry properties for the density correlation 
function (|^5|). In the framework of density functional 
theory the crystal is considered as an extremely inhomo- 
geneous distribution of matter. 

To determine the inverse density correlation matrix we 
use the integral version of the Ornstein Zernike relation 



'5(ri -rs) = I d% {5p{v^)5p{v^)) C{v^,v^) , 
C ra.ra) = c r2,r3 . 



(48) 



The first term of the inverse C(r2,r3) is the ideal gas 
contribution, whereas the second part, the direct cor- 
relation function c(r2,r3), is the contribution from the 
excess free energy, and describes the interactions. More 
precisely, c[r2, ra; n(r)] is a functional of the equilibrium 
density, and is obtained by the second functional deriva- 
tive of the excess free energy J-'^^ with respect to density 
n(r) 

c[ri, ra; n(r)] = /j/Q^ = c\v^^ n; n(r)] . (49) 

In the following manipulations the symmetry expressed 
in (|45|) is used to derive an expression for the inverse den- 
sity correlation function Jgg'(q) in terms of the direct 
correlation function c(ri,r2) starting with the Ornstein 
Zernike (OZ) equation. The left hand side of (|15)) be- 
comes 

y d-Viy d'V3e*(s+i) '-ie^'(s"+'i) '"^ (5(ri3) = l^^gg" • 

(50) 

The right hand side is 

d\x^ d'V3e'(s-Hq)rig-«(g"-Hq)r3x 

X ((5/9(ri)(5p(r2)) C(r2,r3) 
= E'^"g (-G/2-g-q) C(-G-g-q,g" + q) , 

G 

(51) 

which, with Eq. (|47cp . yields 

=Y.^h*{^ + q)'5p(g' + q))-^C(-g' - q, g" + q). 

(52) 



g' 



So the inverse density correlation function is a special 
kind of Fourier transformation of essentially the direct 



correlation function. 

Jgg'(q) = ^C^(-g-q,g' + q) 



y J J \n{ri) 



(53) 



With the definition of the direct correlation function (|49|) 
and its symmetry under interchange of ri o r2 the her- 
mitian property of Jgg'(q) follows 

Jgg'(q) = ^-g.,-g(-q) = (^g'(q))*^ = • (^4) 



B. Invariance under Global Transformations 

One of the fundamental results of density functional 
theorj*22r^ is that the external potential V^^*{r) is a 
functional of the equilibrium density n{r). That is, 
for a given equilibrium density the external potential 
yextfj.-^ = l/'^^*[r; n(f )] is uniquely determined. 

A functional Taylor expansion thus yields 

(5F'=^*(ri) 



(55) 



As the internal state of a crystal is unaffected by a global 
translation and rotation, we now consider the effects of 
such transformations explicitly. This derivation can also 
be considered an invariance principle^i^i, and yields re- 
lations of the direct correlation function for a crystal. 
In the classical approach to elasticity, particle interac- 
tion is described via a potential^^. In that approach the 
consequences of invariance are the conditions for the mi- 
croscopic expressions of the derivatives of the particle 
potential, which ensure the macroscopic Voigt symmetry 
of the elastic constants. We derive analogous results now 
using the direct correlation function. 



1. Translational Invariance 

In the case of a simple translation the transformation is 
given by r' = r -I- s. and the functional Taylor expansion 
yields 

V-^Xv^ + s) V'-Xv^)+jd% ^^(^ ("(^2+8) - n(r2)) . 

(56) 

For an infinitesimal translation — )■ and with a further 
relation from density functional theory 

Ay«^*(ri) 



c(ri,r2) 



'5(ri2) 



(57) 



(5n(r2) ^(-j.^^) 

one obtains the Lovett, Mou, Buf^, Wertheim^ equa- 
tion (LMBW) 



vW[ln7i(ri) + /3y^-*(ri)l = / d\2c{v^,T^)V(;^^n{Y^ 



(58) 
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Without external potential the trivial solution is the one 
with a homogeneous equilibrium density, i.e. a fluid. As 
we arc interested in periodic equilibrium densities, the 
limit of vanishing external field is taken which leads to 
a non-trivial solution. The right hand side is interpreted 
as an effective force on a particle due to interactions with 
the other particles^^. 



Constants of Elasticity 

In this paragraph the dependence on wave vector q 
of Aq^ (q) in the hydrodynamic limit is derived with the 
help of the LMBW equation ([58]) . To do so three con- 
stants of elasticity arc introduced 



Aa^}(q) = ^(5' + q)anl,J*,g{q)ng{g + q)f} 
g'.g 

= (^a/3(q) - iqafJ-isin) + iqpfiaiq) + Qavin) 



(59) 

according to the explicit powers in q. The matrix Ac(/j(q) 
is hermitian, which is a consequence of (j54p . 

We now discuss the three constants of elasticity sep- 
arately, and start with the simplest case, which is the 
term proportional to qaQp- The realness of the equilib- 

n*(r) 



rium density n(r) = J2g n^e'^^'^ = ^ 
yields 



^g g 



g.g 

fcsT 



V 



d\i j d'V2n(ri)n(r2)e-*'i ' 



^(ri2) 
n(ri) 



c(ri,r2) 



where the homogeneous constant v equals 



(60a) 



(60b) 
(60c) 



V 



■j d^-i j d% (?i(ri)5(ri2) - 7i(ri)c(ri, r2)n(r2) 

(60d) 



It also follows, that 



(60e) 



is even in q. The fact that G R and that it has only 
even powers in an expansion in g is a consequence of the 
ri <-> Y-2 symmetry. One further interesting fact is, that 
the equation for v reduces to the inverse compressibility 
K of a fluid'^- for n(r) = and c(ri, r2) = c(ri2) 



V 



[N-nlV j c(ri2) A 12 ) . ((4.27) in (32|) 



The next term ^a{f\} is manipulated with the help 
of the gradient of the equilibrium density VQ,n(r) = 



X]g Wa'Tige^^''^ , and the LMBW equation in the limit of 
vanishing external potential 



A*a(q) 



n*,J*, Jqjrigiga 



'"g'^g'g 

g.g' 



(61a) 



V 



J d^-i j d'V2n(ri)Van(r2)e-*i "-^^ 



'^(ri2) 
?i(ri) 



- c(ri,r2) 



(61b) 



V 



d\i j <i'^r2n(ri)Van(r2)c(ri,r2) 

(61c) 



-iq ri2 



iiJ-apqp + C('?^) , 



(61d) 



where the second rank tensor describing the long 
wavelength limit equals 



^""^ < d\ / A2n(ri)V„n(r2)c(ri,r2)ri2„ 



V 



(61e) 



For a crystal with inversion symmetry it can be shown 
that the correction in the expansion of /iQ(q) is 0{q^). 
The realness of the gradient of the equilibrium den- 
sity Van(r) G R, i.e. V an{T) = X]g «5a"-ge*^''' = 
X)g ~i9an*^P~^^'^ 1 together with the LMBW equation is 
used for the last term Xapiq) 



^a/3(q) = ^ -igLn^'J^'gicihgigfi 

S:S' 

knT 



(62a) 



g.g' 

' / d'V2V„7i(ri)V^n(r2)c(ri,r2) 



(62b) 
(62c) 



V 

1 - e-'i '"!^^ 

~ Kp-iSq-^qs + 0(9**) 
where the fourth rank tensor equals 

Ki3-tS=^^^jd\jd\'^an{Yx)Vpn{Y:i)c{Tx,r-2)ri2^^ri2^s- 

(62d) 

Obviously, one also finds 

KM) = ^}M)- (62e) 

Again it can be shown that due to the ri o symme- 
try, the expansion in q has only even powers, and that 
AQ/j(q) G R, so = A/3Q(q). 

Note also, that for the expansion to be valid, the direct 
correlation function c(ri,r2) is assumed to be of short 
range in the difference vector ri — y-2. 

To sum it up, it was shown in this paragraph that 
limq_>.o AQ^(q) is indeed second order in q, and this was 
derived with the LMBW equation, which is a conse- 
quence of translational invariance. 



AQ/3(q) ~ Aa^^a^^g^ + qafJ-p-yq^ + qisfJ-ajq-y + la^qp 



(63) 
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The second term in Sec. IIII Bl with undetermined q 
dependence was 



\4(q) =t7„t J*— = V(g' + (7)„n;,J*, (q)r 



no 



no no ^ 



no 



(64) 
(65) 



Thus, the momentum equation (j36]) indeed contains a 
term proportional to the gradient of the defect density. 



2. Rotational Invariance 

As translational invariance was the reason be- 
hind the correct q dependence of hmq_>o (q) and 
hmq^oVQ(q), the consequence of rotational invariance 
is now considered. It will be shown, that it yields sym- 
metries in the indices of the constants of elasticity fiafj 
and Xa/s-rS- 

An infinitesimal rotation is given by 

r' = r + de xr + 0{5e^) . (66) 
Thus the first order term of the expansion in S9 is 

^2) 



d% , / '^ J6)xr2-V(^)n(r2). 
dn{r2) 

(67) 



With invariance of the scalar triple product under cyclic 
permutations and an arbitrary S6, one finally ends up 
with a rotational analog of the LMBW equation^ 



rix 



v(i) 



lnn(ri)+/3F"''*(ri) 

= y"d''r2c(ri,r2)(r2xV(2)n(r2)). (68) 



This equation may be interpreted as a balance of effective 
torques in analogy to the balance of forces. 



Symmetry of Constants of Elasticity 

The results for translational, Eq. ([55]). and rotational, 
Eq. invariance can be combined to understand the 

index symmetries of the constants of elasticity. The dif- 
ference (ri X ([58]) ) — ([68]) . which is vahd for any F'^^*(ri), 
yields 

yd''r2c(ri, r2)ri2,„V^''n(r2) =|d'^2c(ri, r2)ri2,;3VPn(r2). 

(69) 



Integrating the last equation with J (i'Vin(ri) leads 
to 







V 



I Ain(ri) J d%c{r^,r2) 



(2). 



(70) 



This is nothing but the statement that fiap is symmetric 
in its indices 



Ma/3 = fJ-Pa 



(71) 



In the same manner, as translational and rotational 
invariance led to a symmetric matrix /Xq.^, the symme- 
try of the indices of A^^-^^ can be addressed. So far it 
is known that Xap-fS = Xpa^yS (consequence of symme- 
try of Aa^(q)) and Xa/s-yS = Xaps-y (symmetric combina- 
tion in 7 (5 in definition (|62p ). Repeating analogous 
arguments concerning the symmetry of Hap, one finds 
that Xap-ys is symmetric under the pairwise interchange 
a/3 o 7(5 (see [2^ for details) 

=> Xap.yS = XjSaP ■ (72) 



V. RESULTS AND DISCUSSION 
A. Summary of the Derived Equations of Motion 

Because the results for the equations of motion are 
spread over different sections, it appears useful to collect 
and list them. 

Starting with the conserved (neglecting for simplicity 
energy conservation) and symmetry restoring fields, we 
showed that the ansatz Eq. (|34aP leads to an exact solu- 
tion of the (for simplicity dissipation-less) Zwanzig-Mori 
equations ([25| . if the seven (d + 1) coarse grained fields 
satisfy the following (because of our use of the fiuctu- 
ation dissipation theorem necessarily) linear equations; 
they hold for all q in the first Brillouin zone. 

Mass density times the time-derivative of the displace- 
ment field equals the momentum density field 



mno dtSuaiq, t) = 5ja{q,t) 



(73a) 



Density fluctuations arise because of the divergence of 
the displacement field and defect density fluctuations 



i5n(q, t) — —inoqa6ua{ci, t) — (5c(q, t) 



(73b) 



Mass is conserved, which connects density and momen- 
tum density fluctuations 



mdtSn{q, t) + iqaSjaiq., t) ~ 



(73c) 



Momentum density, displacement and defect density field 
arc coupled in 

dtSja{(l,t)^-Aai3{q) Sup{q, t) - Va(q) Sc{q,t) 

= - (Aa^i(q) -«gQAi/3(q)+ig,3Ma(q)+gQi^(q)9,3)<5M^i(q, i) 

- — (Ma(q) - iqai^{q))Sc{q,t) . (73d) 
no V / 
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Also the wave equation for the momentum field holds 
mnQdl5ja{(l,t) = -Aa/3(q) 5jp{ci,t) 

= (^ai3ici)-iqafJ-pici)+iqpfJ-*ai(l) + Qav{q)qp^dji3{q,t) . 

(73e) 

In order to recover the momentum wave equation 
Eq. (j73ep by taking a derivative with respect time of 
Eq. (|73d|) and using Eq. (|73ap . the defect density has 
to be constant 



5t(5c(q,t) =0 



(73f) 



Taking the time derivative of Eq. (|73ap and using 
Eq. (j73dp . one sees that the defect density plays the role 
of an inhomogeneity in the wave equation of the dis- 
placement field, which otherwise contains the identical 
constants of elasticity as the momentum one. 

In the hydrodynamic limit, where g ^ 0, the elastic 
coefficients in Eqs. (|73d|) reduce to the following expres- 
sions 

iy{q) = u + 0[q^) , (74a) 
/ia(q) = iiiafsqfi + C(g^) , (74b) 

AQ/3(q) = Kpjsqjqs + 0(9") , (74c) 

with the following symmetries 



Ma;3 = M;3a , 

For later comparison this summary is completed by giv- 
ing the momentum equation in the hydrodynamic limit 

dtSj aiq, t) = iqp (lySap + ^lap) 

V /no 

- Wai3^5 qpqs Su^{q, t) , (76) 
with a wave propagation matrix Wap-yS 

(77) 

which is symmetric in a o 7 and aj /36 according to 
Eqs. (|75p . and symmetrized in /? -f-s- 5, as both indices are 
summed over. 



B. Phenomenological Theory 

For the sake of easy comparison it appears worthwhile 
to summarize the results from thermodynamics and clas- 
sical elasticity theory in order to compare with our mi- 
croscopic expressions. Especially of interest is to ver- 
ify that our results obey the symmetry relations derived 



within the phenomenological approaches. The deriva- 
tion of elasticity theory within nonequilibrium thermo- 
dynamics considering the conserved densities (mass, mo- 
mentum, and energy) and the broken symmetry variable 
displacement field can be found in the literature^ii and 
the result will just be quoted for the reversible, isother- 
mal, and linearized case. 

Let Uai3 = \ (VtjM^ -I- V jjUa) bc the (symmetrized) 
gradient tensor of the displacement field, which agrees 
with the strain field in the considered small deformation 
limit. The first law for the free energy density (per vol- 
ume) as functions of density n and strain Uap is 



df = iidn + hapduaf) , 



(78) 



with chemical potential ^ = jf- 



and hn 



df 



dUa 



the stress tensor at constant density. Note that we keep 
the temperature T constant throughout. 

For the linearized equations of classical elasticity the- 
ory, one requires the isothermal free energy as an expan- 
sion around the equilibrium value hq, which is given by 



/ = /o + 



df 



df 



dUaf 



5n^ 



d^f , 

-2—— On-Ua/s- 

anouap 



duapdujs 



(75a) = /o + fi°Sn + h'l^pUap 



-A 



5n 



no 



(79) 



(80) 



With the equilibrium values of the chemical potential /x*' 
and of the stress tensor at constant density h^^. The 



thermodynamic derivatives are: A 



no 



13- 

an inverse 



ndUaji 



a matrix of coupling 
the elastic coefficients 



compressibility. Bap 

constants, and CTo^a — a 5 — 

at constant density. Due to rotational invariance, i.e. a 
symmetric strain field Uap, the thermodynamic deriva- 
tives obey certain symmetries: the matrix Bap = Bpa 
is symmetric, with up to six independent coefficients de- 
pending on crystal symmetry; and the elastic constants 
(additional symmetry due to definition as second deriva- 
tives) have the Voigt symmetry with a maximum of 21 
independent elastic coefficients. 

The equations of motion derived from microscopic 
starting point in the previous sections contain the defect 
density as fiuctuating variable. Thus it is convenient to 
introduce the defect density as thermodynamic variable 
using the connection between particle and defect density, 
Eq. dSSD, 



5n 



-no Uac 



5c 



(81) 



Changing thermodynamic variables from density to de- 
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feet density gives for the free energy 

= —^dc + (Ja/sdUap , 



with fJ,{c,Uai3) = 

stant defect density 



(82) 
(83) 

and the stress tensor at con- 



"q/3 



Of 



dUaf: 



hal3 - noflSaf}^ 



(84) 



Based on the above thermodynamic expressions, the 
phenomenological equations of motion for a crystal can 
be presented 



dtmSn{c[,t) = -iqaSja{(l,t) 
1 



dtSua{q,t) 
dtSja{ci,t) 



mriQ 



-<5ja(q,i) , 



(85a) 
(85b) 
(85c) 



which express mass conservation, that the time derivative 
of the displacement is the momentum density divided 
by the average mass density, and that the (conserved) 
momentum density changes because of stresses. 

To obtain the desired linear equations of elasticity the- 
ory for this set of variables starting from Eq. (|85|) . the 
partial derivatives of the stress tensor with respect 
to the defect density c and Uap are required. Straightfor- 
ward differentiation and use of the expansion of the free 
energy gives 



df 


df 




dn 




c dUal3 


n dn 


Ua,f3 dUaP 



Sn 



A 5n 
no no 



no 



(86) 



Now everything is in place to state the phenomenolog- 
ical equations of linearized, isothermal and dissipation- 
less elasticity theory with which to compare our micro- 
scopic results. With the change from 6n to Sc the hydro- 
dynamic equations of motion are 



dtSc{q,t) 
dt6ua{q,t) 



0, 



mno 



-6ja{q,t) , 



(87a) 
(87b) 



dtSja{q_,t) = i{A5ap - Bap)qp 



(5c(q,t) 



no 



-QpqsiCap-ys- Bapdjs- SapBjs+ ASap6jsj 6u^{q, t) . 

(87c) 

In the last equation the elastic constants at constant de- 
fect density appear 



Cap^S 



d'f 



dUapdUyS c 

^aPjS ~ BapSjS — 5apB^s + A5ap5^s , (88) 



and a combination of thermodynamic derivatives A5ap — 
Bap, which are based on a Maxwell relation. A re- 
lated combination also showed up in the coefficient in 
the elasticity equations of motion connecting the prefac- 
tors of the displacement and of the defect density field, 
see Eq.(l76l). 

Comparing the classical equations of elasticity theory 
with Eq. (j76p . derived from the Zwanzig-Mori equations 
in the hydrodynamic limit, we can conclude complete 
agreement considering the wave vector dependence, but 
the issue of identifying the microscopic expressions with 
the elastic constants remains open. 



C. Identification of Elastic Constants 

So far wc have shown, that AQ,^(q) cx and the sym- 
metries of the constants of elasticity Xap-yS and ^ap ■ The 
piece, which is still missing, is how these constants are 
related with "the elastic constants" Cap^s- 

As a first observation, the term in front of Jc(q, t) in 
Eqs. ([75)) and (|87cp implies that the coefficient i/ equals 
(up to a constant Co) the thermodynamic derivative 
which was abbreviated as A 



2 d\f 



Co . 



(89) 



Also, the coupling of the density and strain fluctuations, 
abbreviated as B^p^ is given by the matrix —flap 



no 



d'f 
dnduap 



-fJ-ap + CoSap ■ 



(90) 



For the term in front of Suy{q,t) in Eqs. ([7S)) and 
(|87cp the indices and 6 are summed over. Conse- 
quently the fourth-rank tensor of wave propagation coef- 
ficients Wap^s, which was defined already symmetrized, 
has to be compared with symmetrized elastic constants 
^{Cap-fS + CaS-fp) = Wap-js- This yiclds the rele- 
vant combination^i^d^ for the elastic constants Cap-yS in 
terms of wave propagation matrices Wap-fS, or, respec- 
tively, constants of elasticity v, flap , and Xap-yS 

Cap^S ~ (WaPjS + WpajS - Wa-yPS^ 

= (^Xa^pS + XpjaS — Xap'yS^ +'^5ap5jS 

+ flapSyS+fJ-ysSaP ■ (91) 

Several interesting results for this combination might be 
noted: the set of three independent Xap-yS, which are not 
related via Voigt symmetry, occur in the combination for 
the elastic constants. The combination of flap and ly are 
only in pairs of the indices a/3 and jS; there is no Voigt 
symmetric term fiaS^p-y + fips^a-y + fJ-a-ySps + fip-ySaS or 
i'{5a'ySps + SaS^P'y)^ which, for an isotropic solid, corre- 
sponds to the combination of the shear modulus. 
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The elastic constant at constant density are thus given 
by the matrix A defined in Eq. (|62p via 



dUafiUyS 



(92) 



The derived results for the elastic constants in terms 
of the direct correlation function parallel other known 
expressions for quantities characterizing broken symme- 
tries in terms of c(ri,r2). The Triezenberg-Zwanzig 
expression^ for the surface tension between gas and liq- 
uid phase of a phase separated simple system contains 
the equivalent quantities as our results, namely the direct 
correlation function and the average density profile. For 
the surface tension, Kirkwood and BufS^ gave another 
equivalent expression in terms of the actual interaction 
potential and the density pair correlation function. For 
the elastic coefficients familiar results in terms of the par- 
ticle interaction potentials can be found in the classical 
textbooks2ii2^, yet only for the case of ideal crystals and in 
the limit of low temperature where particles fluctuate lit- 
tle around the lattice positions. Shortly, the connection 
can be established via the symmetrized wave propagation 
coefficients 



with ^ajihj) = Va V^^^$(R', Rj), which contains the 
actual potential ^> and is evaluated at the equilibrium po- 
sitions R. Interestingly, we can recover Eq. (^5)) using the 
mean-spherical approximation c(ri, r2) = — /3$(ri,r2) 
and appropriately coarse-graining^^. However, we are 
not aware of results equivalent to ours containing the ac- 
tual potential and the pair correlations functions at finite 
temperature in non-ideal crystals. 



D. Displacement and Defect Density Fields 

After the identification of the coefficients appearing 
in elasticity theory, it is worthwhile to turn to the mi- 
croscopic definition of the displacement field which re- 
sulted from following the standard approach to Zwanzig- 
Mori equations of broken-symmetry systems. For crys- 
tals of cubic symmetry, where Afap = ■N' 6ap simplifies in 
Eq.dH]), it is 



g 



(94) 



Importantly, this relation allows to determine the dis- 
placement field purely from measuring density fluctua- 
tions. No reference lattice is required. Thus, this formula 
can be applied in non-ideal crystals containing arbitrary 
concentrations of point defects like vacancies and inter- 
stitials. For our equilibrium considerations to apply, the 



point defects should be mobile and diffuse during the 
measurement, even though defect diffusion is neglected 
yet in our dissipation-less formulation. Equation (|4T|) 
was first ingeniously formulated by Szamel and Ernst^, 
who also started with the Fourier components of the den- 
sity, but then changed to the usual set of hydrodynamic 
variables. Due to their change of hydrodynamic variables 
their results differ from ours in the interpretation of i> as 
inverse isothermal compressibility, and in the neglect of 
the coupling term (Su'^Sn); see the Appendix. Because of 
this, Szamel in a continuation paper— concluded incon- 
sistencies to phenomenological elasticity theorj*i, which 
are absent in our results. 

Also the result for the density fiuctuation appearing in 
the equations of elasticity theory is noteworthy 



(95) 



This coarse-grained density fiuctuation differs from the 
microscopic density ffuctuation defined in Eq. ()14|) . One 
of the consequences is that the correlation function of the 
coarse-grained density is not simply related to the gener- 
alized structure factor iS'g(q) defined in terms of micro- 
scopic density fluctuations in Eq. (jTZ]). In the limit of 
zero wave vector, the structure factor at g = reduces 



(93) to the isothermal compressibility, 5'g=o(q) 



kfiTu 



while the coarse-grained density ffuctuation function re- 
duces to the thermodynamic derivate {Sn*{q)Sn{q)) 

'■^^ = VkuT (0)"' (for Co = 0). 



VI. CONCLUSIONS AND OUTLOOK 

The definition of a displacement field is central to the 
description of crystal dynamics. Yet, for non-ideal crys- 
tals containing point defects, it had been lacking. We 
provide the first systematic derivation of a microscopic 
expression for the displacement field in terms of den- 
sity fiuctuations with wavelengths close to the recipro- 
cal lattice vectors. We also find that the coarse grained 
density field of elasticity theory differs from the (naively 
expected) small wavevector limit of the microscopic den- 
sity. These expressions lead to microscopic formulae for 
the elastic constants of a crystal given in terms of the 
direct correlation functions. A discussion of the symme- 
tries of the direct correlation functions recovers the (re- 
quired) symmetries of the elastic coefficients for general 
crystals. Complete agreement with the phenomenolog- 
ical description given by the linearized, isothermal and 
dissipationless elasticity theory is achieved. 

A generalization of the approach to include energy 
fiuctuations and dissipation is possible; as are exten- 
sions to other broken symmetry systems, like quasicrys- 
tals and liquid crystals. Closure approximations for 
the direct correlation functions^ will enable quantita- 
tive evaluation of the derived formulae. A generalization 
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of the theory is required for the inclusion of topological 
defects^, which destroy the order parameter rig. This 
could then be compared with the continuum theory of lat- 
tice defects^i^. Colloidal crystals would provide model 
systems^"— where the theory can be tested. 
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Appendix A: Conventional Set of Variables 

For comparison wc outline in this appendix how the 
Zwanzig-Mori equations of motion arc derived with the 
conventional set of variables, i.e. the density (5n(q, t), the 
momentum density 6ja (q, t) , and, as broken symmetry 
variable, the displacement field dua{q_,t). In order to 
explicitly calculate correlations containing density and 
displacement field, we use their expressions in terms of 
the microscopic density fields given in Eqs. (PT|) and (|^^ . 
The main result of this appendix is the interpretation of 
the terms AQ;3(q), /^a(q), and J^(q). 

The Zwanzig-Mori equations of motion Eq. with 
the slow variables (SAi) G {Sn, 5ja, ^Uq} contain the non- 
vanishing matrix elements of the Liouville operator 

(5i;(q)£5n(q)) = -qc^noVkeT , (Al) 
{53l{ci)C5up{ci)) = ~i5a.i}VkBT . (A2) 

While the former macroscopically follows immediately 
with Eq. (|73cp from the equipartition theorem, Eq. ([T5| , 
the latter additionally requires identifying the time 
derivative of the displacement as momentum density (di- 
vided by mass density), Eq. (|73ap . Under the assump- 
tion, that the microscopic fluctuations Sp{g + q) may be 
replaced by the hydrodynamic ones (5ng(q), 

Sti( ) 

6p{g + q) ~ Srigiq) = -iUggaUaiq) + rig , (A3) 

no 



both matrix elements arc rcderived from Eq. (jlSp . In the 
first case this leads to 

(A4) 

and in the second case to 

{^jai^)^Sup{q)) = -iM^^^n*^g^{ga + qa)n^VkBT . 

(A5) 

Due to the definitions of TV^^ in Eq. (jH]), A/'o in Eq. 
and the orthogonality Eq. (j40[) . both summations red- 
erive Eqs. (|XT|) and (|X2|) . 



Because of time reversal symmetry the Zwanzig-Mori 
equations of motion Eq. contain the following non- 
vanishing isothermal correlations: the equipartition the- 
orem for the momentum density, Eq. (jl6p : the correla- 
tion of the coarse-grained density, {6n* {q)Sn{q)); the dis- 
placement correlation function, (Jri* (q)(5M/3(q)), which is 
the inverse of the dynamical matrix at constant density; 
and a coupling between the displacement field and the 
coarse-grained density fluctuation, (i5m* (q)5n(q)). The 
latter three correlations are calculated most easily to- 
gether, as follows with the assumption Eq. (|A3[) . We 
now consider the consequences for the two-point correla- 
tion functions 

(V(g + q)W + q)) = (A6) 

^Sn*{q)6n{q) , 



igiig. 



(5a.9^« (q)"/3(q)) + (- 



-) 



+ Wc^Mq) ) - «( up{q))g^ 

no no 

This is then inserted in the Fourier transformation of 
the Ornstein-Zernike equation ([M)) where the definition 
of the inverse density correlation (|53)) enters. With the 
same manipulations as in chapter llVB H one obtains 

„ „ , N,<5n*(q) , „ , , ,5n*(q)(5n(q) , 
VkBT=-p^{q){ ^^^(q))+z.(q)( ^^^ ), 



no 



(A7a) 



= /i/3(q)«(q)"/3(q)) - Kq)«(q)^), 



no 



(A7b) 

^ , .,Sn*(q) , .. ^, .,Sn*{q)Sn{q)^ 
0^-Xpsiq){^^U5iq))+IJ,*p{q}{ ' ^ 

(A7c) 



no 



vkBTSafi^Xp5{q){Kiqhsiq)) - /ij(q)«(q; 



Sn{q) 



no 



(A7d) 



From this set of equations, the three desired isothermal 
correlations can be read off. The correlation of the con- 
served density or coarse-grained structure factor is given 
by the inverse of i'(q) 



(Sn (q)(5n(q)) = ^^^^ 



(A8a) 



The displacement correlation function (at constant den- 
sity) is given by the inverse of the matrix Xap 



{Su*M)Sup{q)) = VkBT{X^0{q))-' 



(A8b) 



which shows that the displacement fluctuations are long- 
rangcd, (q)(5w^(q)) cx for g — >• 0. Lastly, the 
coupling between the displacement field and the coarse- 
grained density fluctuation is given by the inverse of 
fJ-aiq) 



{Su*M)Sn{q)) = 



-npYkBT 

Ma(q) 



(A8c) 
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Note that Szamelii assumes that this last correlation is The result is, with Eq. (|35|) . equivalent to Eqs. ([73)) . 
0(q), i.e. negligible in the hydrodynamic limit, while we 
find that it grows like l/q and can not be neglected. 

Finally the Zwanzig-Mori equations of motion of a 
crystal in terms of the coarse-grained density, the mo- 
mentum density, and the displacement field are 

dtSn{(i,t) = ~i—qaSja{q,t) , (A9a) 
m 

dt5ja{(l,t) = ~iqa'^^Sn{q_,t) ~ Ac<^(q)5u^(q, t) 
no 

+ iqc.^^|3{^)Sufi{q,t) + ^^^^^Sn{q,t) , 
no 

(A9b) 

dt6ua{q,t) —^6ja{(l,t) . (A9c) 
mno 
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